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ABSTRACT 

Brady rhizobium elkanii is successfully used in the formulation of commercial inoculants and, together with 
B. japonicum, it fully supplies the plant nitrogen demands. Despite the similarity between B. japonicum 
and B. elkanii species, several works demonstrated genetic and physiological differences between them. In 
this work Representational Difference Analysis (RDA) was used for genomic comparison between B. 
elkanii SEMIA 587, a crop inoculant strain, and B. japonicum USDA 1 10, a reference strain. Two hundred 
sequences were obtained. From these, 46 sequences belonged exclusively to the genome of B. elkanii 
strain, and 154 showed similarity to sequences from B. japonicum genome. From the 46 sequences with no 
similarity to sequences from B. japonicum, 39 showed no similarity to sequences in public databases and 
seven showed similarity to sequences of genes coding for known proteins. These seven sequences were 
divided in three groups: similar to sequences from other Brady rhizobium strains, similar to sequences from 
other nitrogen-fixing bacteria, and similar to sequences from non nitrogen-fixing bacteria. These new 
sequences could be used as DNA markers in order to investigate the rates of genetic material gain and loss 
in natural Brady rhizobium strains. 

Key words: Brady rhizobium; Genome comparison; Nitrogen fixation; Representational difference 
analysis. 



INTRODUCTION 

Soybean [Glycine max (L.) Merrill] can establish effective 
biological nitrogen fixation (BNF) symbiosis with species of 
fast-growing rhizobia (13, 23, 25) as well with species of slow 
growing bradyrhizobia, like Brady rhizobium japonicum (27), 
B. elkanii (30), and B. liaoningense (50). In theory, the 
nitrogen sources available for soybean cultures are nitrogen 
fertilizers and BNF (23). In Brazil, for example, the use of 
nitrogen as fertilizer is not recommended for soybean, as the 
inoculation with nitrogen-fixing bacteria fully supplies the 
plant nitrogen demands. Besides the saving generated annually 
by the inoculation - about US$ 3 billion per year (24), the 



practice avoids environmental pollution, very common when 
nitrogen fertilizers are employed. Since soybean Brazilian's 
fields lacked native bacteria able to nodulate this leguminous, 
the production of inoculants in this country started in the 1960s 
with North American strains (12). Nowadays, only four B. 
japonicum and B. elkanii strains have been used in the 
formulation of Brazilian's commercial inoculants and they 
have resulted in an established population in most soils 
cropped with soybean (15). Of those recommended strains, 
SEMIA 587 and SEMIA 5019 belonged to the B. elkanii 
species (42), while SEMIA 5080 and SEMIA 5079 belonged to 
the B. japonicum species (32). 

In the 1980s the genus Brady rhizobium was originated 
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from the slow- growing specie Rhizobium japonicum (27). In 
the beginning the soybean-nodulating B. japonicum was the 
only specie described. To date, five additional species have 
been vaUdly named in this genus, and besides B. japonicum 
other two species are able to nodulate Glycine max, B. elkanii 
and B. liaoningense (49). Despite the similarity between B. 
japonicum and B. elkanii, years of research demonstrated 
genetic and physiological differences among them. Boddey and 
Hungria (7) when studying 40 soybean Bradyrhizobium strains 
showed remarkable differences between these two species 
according to the resistance to antibiotics, where B. elkanii 
SEMIA 587 was resistant to kanamycin, rifampicin, 
spectinomycin, streptomycin, carbenicillin, chloramphenicol, 
nalidixic acid, tetracycline and erythromycin, while B. 
japonicum USDA 110 was resistant only to low levels of 
rifampicin and spectinomycin. These authors also found 
differences concerning to the production of indole acetic acid 
(lAA): B. japonicum strains accumulated between 4.88 to 7.08 
|iM of lAA ml"\ while B. elkanii strains were able to 
accumulate about six times more (44.36 [iM of lAA ml"^). 
Genetics studies like the sequencing of nifDK and nijE genes 
(34) and the 16S rRNA gene (51), together with hybridization 
analysis with the hup gene (35), nifDH genes, and homologous 
sequences of nodD gene (44) also helped in the final division 
of B. japonicum and B. elkanii made by KuykendaU and 
coworkers (30). However, up to now the majority of published 
works about the genus Bradyrhizobium are related to B. 
japonicum strains. Nowadays, B. japonicum USDA 110 (28) 
and Bradyrhizobium sp. BTAil (17) are the only 
Bradyrhizobium species to have the whole genome sequenced. 
Concerning B. elkanii species little is published and known 
about it, especially in the field of genetics and genomics. 
Basically nif, nod, recA, glnll, and ribosomal genes sequences 
are found in the public databases. 

The original Representational Difference Analysis (RDA) 
protocol was developed for detection of genomic differences 
between two genomes (33). The methodology has been 
modified and employed as a sensitive method for cloning of 
differentially expressed genes (1, 4, 8, 31). Microbiologists 



also succeeded to reveal genomic differences by comparing 
bacterial strains belonging to the same species with contrasting 
properties (6, 18). 

Due to the ecological and enviroimiental role of the BNF, 
the economy and efficiency of inoculants and their massive 
usage over soybean crops, it is urgent a better understanding of 
the B. elkanii genetics. To obtain such understanding B. elkanii 
SEMIA 587 and the reference strain B. japonicum USDA 110 
were submitted to RDA experiment. 

MATERIALS AND METHODS 

Bacterial strains, growth conditions and DNA extraction 

The rhizobial strains used in this study were the reference 
strain Bradyrhizobium japonicum USDA 110, and the crop 
inoculant strain Bradyrhizobium elkanii SEMIA 587; both 
strains were supplied by Funda9ao Estadual de Pesquisa 
Agropecuaria (FEPAGRO, RS, Brazil). Each strain was grown 
on Yeast Mannitol Broth (10 g mannitol, 0.5 g potassium 
phosphate, 0.2 g magnesium sulphate, 0.1 g sodium chloride, 
0.5 g yeast extract, water to one Uter, pH 6.8) for seven days at 
28°C in an orbital shaker. Genomic DNA from bacteria was 
extracted according to basic molecular biology protocols (43). 

RDA 

The technique was modified from that described by 
Tinsley and Nassif (46). DNA from B. japonicum strain USDA 
110 was used as driver DNA and DNA from B. elkanii strain 
SEMIA 587 was used as tester DNA. Fifteen |ig of driver DNA 
were used for the whole RDA experiment. One ng of tester 
gDNA was completely digested with 5'aM3AI, resulting in a 
population of fragments with 200-300 pb, that were purified by 
colunm (Wizard® SV Gel and PGR Clean-Up System - 
Promega) and ligated to 500 pmol of each adapters NBam 12 
(GATCCTCCCTCG) and NBam 24 (AGGCAACTGTGCTAT 
CCGAGGGAG). The Oscilating Phenol Emulsion 
Reassociation Technique [OsPERT, (10)] was used for the 
subtractive hybridization. In the first round 50 ng of tester 
DNA were mixed with 5 |ig of driver DNA, resulting in a rate 
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of 1:100, and held for 15 min at room temperature to denature 
in a final volume of 34.6 |j1 of 0.3 N NaOH. Solution was 
neutralized with addition of 5.4 pi of 3 M MOPS, 50 |j1 of 1.5 
M NaCl, 2 mM EDTA. Ten pi of buffer-saturated phenol pH 
8.0 was added and mixed. The sample was heated at 65°C for 1 
minute and put on ice to generate a phenol emulsion. Then it 
was placed in a thermal cycler (PGR Express ThermoHybaid 
Thermal Gycler) programmed for 99 cycles of 15 seconds at 25 
°G and 15 seconds at 65 "G. The phenol was removed with 
chloroform and the supernatant was purified by column. Five 
|j1 of the hybridization was added to 100 pi of reaction mix 
without oligonucleotides (IX buffer, 0.25 mM of each dNTP, 
25 mM of MgGl2, 2 U Taq Polymerase - Invitrogen) and 
incubated for 3 min at 72°G to fill in the ends of the reannealed 
fragments. After denaturation at 94 °G for 5 min and addition 
of the oligonucleotide NBam 24 (100 pmol) the hybridized 
fragments were amplified by PGR (30 cycles of 1 min at 94 °G, 
1 min at 62 °G, and 1 min at 72 °G followed by 5 min at 72 °G 
- PGR Express ThermoHybaid Thermal Gycler). The PGR 
products were again purified by column. An aliquot of these 
purified fragments was separated for a 1:10 dilution. The rest 
of them were treated with 1 U of Mung Bean Nuclease 
(Promega) for 15 min and column purified. Two PGR were 
performed, one reaction used 2 pi of the 1:10 dilution as 
template and the other reaction used 10 pi of the Mung Bean 
treated fragments as template. These PGR were performed as 
above described without the filling step. Next, the adapters of 
the differential products 1 (DPI) were removed by Sau3Al 
cleavage and a new set of adapters was ligated to DPI [JBam 
12 (GATGGGTTGATG) and JBam 24 (AGGGAGGTGGAGT 
ATGGATGAAGG)]. The new adapter hgated DPI (5 ng) was 
used in the second round of subtractive hybridization with 5 pg 
of driver DNA in a rate of 1:1000, and selective amplification 
was performed as in the first round. The third and final round 
was performed as described for the previous rounds but with 
adapters RBam 12 (GATGGTGGGTGA) and RBam 24 
(AGGAGTGTGGAGGGTGTGTGAGGGAG); 500 pg of DP2 
was hybridized with 5 pg of driver DNA resulting in a rate of 
1:10000. 



Construction of subtractive libraries and sequencing 

An aliquot of the final RDA products (DP3 fragments) 
was filled in with Klenow fragment and phosphorylated with 
T4 polynucleotide kinase (Promega). Upon heat inactivation, 
the samples were column purified and hgated to 5mal-digested 
and dephosphorylated pUG18 vector. Escherichia coli XLl 
Blue competent cells were heat-shock transformed with the 
ligation products. White Escherichia coli colonies were 
selected and extracted plasmids were sequenced in both 
directions in the AGT Gene Laboratory (Gentro de 
Biotecnologia, UFRGS, Brazil) using the automatic sequencer 
ABl-PRISM 3100 Genetic Analyzer and Ml 3 universal 
primers. Labeling reactions were performed in a Gene Amp 
PGR System 9700 (Applied Biosystems) thermocycler. The 
sequences were analyzed with the Staden Package software (v. 
1.7.0, http://www.sourceforge.net) . The final resulting 
sequences were compared to the GenBank database using the 
BLASTN program. 

RESULTS 

RDA procedure 

Three rounds of subtractive hybridization and kinetic 
enrichment were performed using DNA from B. elkanii 
SEMIA 587 as tester and DNA from B. japonicum USDAUO 
as driver. As expected, exponential amplification of fragments 
specific to the tester DNA was obtained (see below). Linear 
amplification of fragments that were present on either DNA 
populations (tester and driver) was also observed. Several 
works had suggested the use of Mung bean nuclease to 
eliminate or reduce the unwanted linear amplification during 
the amplification step (38, 45) while others suggested that a 
1:10 dilution is enough for the achievement of this proposal 
(2). Both alternatives were used in the present work during the 
PGR step. As similar amplifications results were achieved 
using either Mung bean or 1:10 dilution as template, both 
amplifications products were mixed and used as one in the 
following RDA rounds. 
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BLAST Results 

Several genetic differences between DNAs from B. 
japonicum USDAllO and B. elkanii SEMIA 587 strains were 
successfully isolated by RDA methodology. After the three 
rounds of subtractive hybridization and kinetic enrichment the 
resultant products were used to construct a RDA library. The 
pUC18 cloned fragments were sequenced, analyzed by Staden 
Package and screened for similarity with those available on the 
GenBank database using BLASTN program. As the complete 
genome sequence of B. japonicum USDAllO is also available 
on GenBank database the RDA library sequences could be 
screened directly against it. 

From the 234 total sequences obtained, 52 were redundant. 



One consensus sequence grouped 1 1 redundant sequences (this 
consensus sequence showed no obvious similarity with any 
sequence available in the pubhc database), five consensus 
sequences grouped three redimdant sequences each, and 13 
consensus sequences grouped two redundant sequences each. 
After grouping the redundancy 200 RDA unique sequences 
were obtained. Forty six sequences among these 200 belonged 
exclusively to the genome of the tester strain B. elkanii SEMIA 
587, and 154 showed similarities to sequences belonged to the 
genome of the driver strain B. japonicum USDAllO. The 154 
sequences with similarity to the genome of B. japonicum 
USDAllO were divided in five groups according to their 
identity (Table 1). 



Table 1. Grouping of 154 sequences of B. elkanii SEMIA 587 with similarity to B. japonicum USDA 1 10 genome. 

Group Number of Percentage of total analyzed Similarity "* 
sequences sequences " 

I 2 1 60 to 69% 

II 9 4.5 70 to 79% 

III 19 9.5 80 to 89% 

IV 68 34 90 to 99% 

V 56 28 100% 

'Percentage based on the 200 consensus sequences. Sequences with no similarity to B. japonicum (46) genome represents 23%. 

Similarity of the clone sequence related to their similar sequence in B. japonicum USDAl 10. 



From the 46 sequences with no similarity to B. japonicum 
USDA 110 genome, 39 showed no similarity with sequences in 
pubUc databases and seven sequences showed similarity with 
sequences of genes that coding for known proteins. These 
seven sequences were divided in three groups: similar to 
sequences from other Bradyrhizobium strains, similar to 
sequences from other nitrogen-fixing bacteria, and similar to 
sequences from non nitrogen-fixing bacteria. The group with 
similarity to sequences from other Bradyrhizobium strains was 
composed by two clones: clone 15 presented similarity to the 
gene encoding a putative toxin secretion ABC transporter fi-om 
Bradyrhizobium sp. ORS278, and clone 129 was similar to the 
gene encoding a putative carboxylase Uke RuBisCO small 
subunit from the photoorganotroph Bradyrhizobium sp. BTAl. 



The group with sinnilarity to sequences from other nitrogen- 
fixing bacteria was composed by three clones: clone il50 
presented similarity to the gene encoding a 4-hydroxybenzoyl- 
CoA reductase alpha-subunit of Mesorhizobium loti, clone il70 
to the gene encoding a conserved hypothetical protein of 
Rhodopseudomonas palustris, and clone ii23 to the gene 
encoding a virulence factor MVIN family protein of 
Xanthobacter autotrophicus Py2. The group with similarity to 
sequences from other non nitrogen-fixing bacteria was also 
composed by two clones: clone 165 presented similarity to the 
gene encoding a peptidase M19 of Sphingopyxis alaskensis 
RB2256, and clone 1157 to the gene encoding a conserved 
hypothetical protein of Nitrobacter winogradsky. Table 2 
shows the complete information about these data. 
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Table 2. List of B. elkanii SEMIA 587 clones with similarity to know genes obtained through BLASTN search. 



Group 


Clone 


Gene" 


Species 


Similarities 

/c 


Expected 
value 


GenBank 

accession 












numbers 


kJ-llllllMJ. O LU UUlti 

fir/1 /i^? rhivnhiiivyi 

strains 


15 
129 


Putative toxin secretion ABC 
transporter 

Putative carboxylase, like 
RuBlsCO, small subunit 


Bradyrhizobium sp. 

ORS278 
Bradyrhizobium sp. 
BTAll 


80 
76 


2e'"' 
2e" 


FT624080 
FT624081 




1150 


4-Hydroxybenzoyl-CoA 


Mesorhizobium loti 


75 


le-'^ 


FJ624082 


Similars to other 




reductase alpha- subunit 










nitrogen-fixing 
bacteria 


1170 


Conserved hypothetical 
protein 


Rhodopse udomonas 

palustris 


70 


5e-°^ 


FJ624083 




1123 


Virulence factor MVIN 


Xanthobacter 


69 


9e-'^ 


FJ624084 




family protein 


autotrophicus Py2 


Similars to other 
non nitrogen- 
fixing bacteria 


165 


Peptidase Ml 9 


Sphingopyxis 
alaskensis RB2256 


76 


3e-'« 


FJ624085 


1157 


Conserved hypothetical 
protein 


Nitrobacter 
winogradsky 


76 


le-^^ 


FJ624086 



' In this case each gene identified means the best hit that was obtained through BLASTN search. 



Sequences belonging exclusively to both B. japonicum and 
B. elkanii strains 

Besides the exclusive B. elkanii SEMIA 587 sequences, 
three sequences of B. elkanii SEMIA 587 that matched only 
sequences from B. japonicum USDA 110 and no other 
organism were found. They were: clone 1184 (107 pb) and 
clone 192 (196 pb) that showed 77% and 94% of similarity, 
respectively, with two genes for different transcriptional 
regulatory proteins, and clone 1176 (91 pb) that showed 74% of 
similarity with an intergenic region (74 bp at 3' side of bsl6114 
and 162 bp at 5' side of blr6113, respectively, a hypothetical 
protein and an unknown protein). 

Sequence acquisition numbers 

GenBank accession numbers of described unique tester 
sequences are: clone 15, FJ624080; clone 129, FJ624081; clone 
1150, FJ624082; clone 1170, FJ624083; clone 1123, FJ624084; 
clone 165, FJ624085, and clone 1157, FJ624086. 

DISCUSSION 

RDA strategy 

There are two basic RDA protocols usually performed to 



find genomic unique sequences: The original protocol 
described by Lisitsyn and coworkers (33) and the modification 
of the RDA procedure described by Tinsley and Nassif (46). In 
the standard RDA approach, chromosomal DNAs of both tester 
and driver DNAs are digested to completion with Bglll and 
complementary oligonucleotides are ligated to them. A 
'representation' of both sets of fragments is generated using 
PGR. Following amplification, oligonucleotides from the tester 
are digested off and the tester is ligated to a new ones. Both 
'representations' are used in the hybridization step. In the 
Tinsley and Nassif (46) procedure tester DNA is digested to 
completion with SauiPd to generate a pool of smaller 
fragments than those generated by BglW digestion. 
Chromosomal DNA from the driver strain is sheared 
mechanically by 60 passes through a 30 gauge needle to 
generate fragments in the 3 kb to 4 kb range. >S'aM3AI digested 
tester DNA is ligated directly to oligonucleotides. Tester and 
sheared driver are mixed and hybridized. Caha and coworkers 
(11) made a comparison of these two protocols while studying 
two strains of Vibrio cholerae. Although the two methods can 
be used, the second was suggested to yield a more complete set 
of different products for analyzing genomic differences. In the 
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present work the Tinsley and Nassif (46) RDA protocol was 
used with some modifications. Instead shearing mechanically 
the driver genome by passes through a 30 gauge needle the 
driver DNA was partially digested with the same SauiP^ 
restriction enzyme (43). Fragments ranging from 1 kb to 10 kb 
were chosen to perform the RDA experiment. 

From the 200 unique sequences, out of 234 sequences that 
were obtained, 46 had no similarity to the sequences from the 
driver DNA, indicating that the RDA library constructed was 
not extinguished and more unique sequences could be obtained 
through more sequencing. The method applied in this work 
yielded a proportion of approximately 20% of unique tester 
fragments. 

The option to compare Bradyrhizobium elkanii SEMIA 
587 and Bradyrhizobium japonicum USDA 110, which has its 
whole genome sequenced, had proved to be a valuable choice. 
Since the complete genome sequence of B. japonicum USDA 
1 10 is available in public databases, a direct comparison of the 
data obtained from tester DNA through sequencing could be 
easily performed, avoiding the time and labor consuming 
Southern blot experiment to the validation of the RDA 
procedure. Additional information were also obtained Uke 
those regarding the tester sequences that have some similarity 
with driver sequences (see Table 1), and the identification of 
sequences that belong exclusively to the genomes of 
Bradyrhizobium elkanii SEMIA 587 and B. japonicum USDA 
110 strains. 

B. elkanii SEMIA 587 unique sequences 

The RDA approach revealed 46 sequences without 
similarity to B. japonicum USDA 110 genome. From these, 39 
showed no sinmlarity with sequences in pubhc databases and 
seven showed similarity with sequences of genes that codify 
proteins, some of them with known function. 

The sequence carried by clone 129 was similar to a 
sequence of the gene that codes for a putative carboxylase 
(small subunit) like Rubisco of Bradyrhizobium sp. BTAil. 
This enzyme catalyzes the biological reduction and 
assimilation of carbon dioxide gas to organic carbon and it is 



the key enzyme responsible for the bulk of organic matter 
found on earth. Its biochemical role is chiefly as the central 
part of the Calvin-Benson-Basham (CBB) reductive pentose 
phosphate pathway where it combines CO2 with ribulose 1,5- 
bisphosphate to form two molecules of 3-phosphoglycerate (5). 
Rubisco forms I and II are directly involved in carbon 
metabohsm and have a well recognized autotrophic CO2 fixing 
function. These enzymes are found in a wide range of chemo, 
organo, and phototrophs, occurring in bacteria, algae, and 
plants. Form III Rubisco has only been found in archaea and 
form IV enzymes have been shown to be functional in the 
methionine salvage pathway in many bacteria. (3). 
Bradyrhizobium japonicum USDA 110 strain, which DNA was 
used as a driver, has an organotroph metabolism; it means no 
photosynthesis is made. The putative Rubisco-like gene found 
in B. elkanii SEMIA 587 genome showed a high level of 
similarity with genes that coding for proteins that contain the 
AHS2 superfamily conserved region, which includes the 
Rubisco class I of Bradyrhizobium sp. BTAil, a facultative 
photoorganomixotroph bacterium (5). As this is the first 
evidence of a Rubisco-like gene in the strain of study the 
necessity to evaluate if this gene is functional in B. elkanii 
SEMIA 587 will be very important. 

The sequence carried by clone 15 was similar to the gene 
that codes for the toxin secretion ABC transporter {hlyB Uke 
protein) of Bradyrhizobium sp. ORS278. ABC (ATP Binding 
Cassette) proteins constitute a specific class of ATPases that 
are involved with a wide variety of transport processes, 
including both inward (import) and outward (export or 
secretion), the latter being conserved from bacteria to humans 
(20, 21). The most studied of these transporters are key 
components of the Type I secretion pathway in Gram-negative 
bacteria. The hemolysin B (hlyB) protein of E. coli is required 
in conjunction with hemolysin D (hlyD) and the extern 
membrane protein TolC for the secretion of hemolysin-A 
(hlAy) protein (19, 48). The Type 1 secretion pathway does not 
only export hlyA but also a great range of substrate like, 
proteases, lipases and different toxins. The substrate can be 
small as 19 kDa or even protein complexes of 800 kDa (22). 
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The systems constitute a complex of three proteins, a 
transporter located in the cytoplasmic membrane, a periplasmic 
membrane fusion protein, and an outer membrane channel 
protein. When the complex is assembled a periplasm crossing 
channel is formed (16). 

Clone il50 contained a DNA fragment which sequence 
had showed similarity to the sequence of the gene that codes 
for Mesorhizobium loti 4-Hydroxybenzoyl-CoA reductase 
alpha-subunit. This enzyme catalyzes an important reaction in 
the anaerobic metabolism of phenolic compounds, i.e. the 
reductive removal of an aromatic hydroxyl group. Anaerobic 
metabolism is based on a reductive chemistry, by which groups 
at the ring are often reductively removed and dearomatization 
of the aromatic ring structures proceeds by reductions. One 
central intermediate is benzoyl-CoA, which is dearomatized by 
reduction. Most natural aromatic substrates are hydroxylated 
compounds that require reductive dehydroxylation of phenolic 
hydroxyl groups if benzoyl-CoA is to be formed (9). 

Homolog to Xanthobacter autotrophicus Py2 virulence 
factor MVIN family protein was first described in Salmonella 
typhimurium. This protein has been shown to be an important 
virulence factor for this organism when infecting mouse (29) 
and it is part of a Sinorhizobium meliloti essential operon (41). 
Moreover, the mviN gene in Rhizobium tropici was found to be 
involved in motility (37) and a mutation in mviN gene in the 
soft rot-causing bacterium Pseudomonas viridiflava was shown 
to influence motihty as well as virulence (14). However, the 
biochemical function of the product encoded by mviN and 
similar genes is not fuUy understood (26). A 69% of similarity 
was found between the sequence present at clone ii23 and the 
gene coding for the Py2 virulence factor MVIN from X. 
autotrophicus. It will be interesting to perform mutagenesis 
studies in B. elkanii in order to verify if this gene is also 
involved in the motihty or other process in this bacterium. 

The sequence carried by clone 165 was similar to the 
sequence of the gene that codes for Sphingopyxis alaskensis 
RB2256 peptidase M19 or renal dipeptidase, which is an 
enzyme that in humans is a membrane bound glycoprotein 
involved in the hydrolysis of dipeptides and is well-known as 



the sole metabohc enzyme for penem and carbapenem P-lactam 
antibiotics in mammals. These antibiotics are presently 
prescribed in chnical practice with the dipeptidyl inhibitor 
cUastatin (36). Peptidases are grouped into clans and families. 
Clans are groups of families for which there is evidence of 
conmion ancestry. Families are grouped by their catalytic type: 
A, aspartic; C, cysteine; G, glutamic acid; M, metallo; S, 
serine; T, threonine; and U, unknown. In metalloproteases a 
divalent cation, usually zinc, activates the water molecule. The 
metal ion is held in place by amino acid ligands, usually three 
in number (40). As already mentioned, B. elkanii SEMIA 587 
was found to be more resistant to antibiotics than B. japonicum 
USDA 110 (7). Bacteria in general has a wide range of 
mechanisms to resist the actions of noxious agents: enzymatic 
reactions like P-lactamase, which hydrolyzes the P-lactam ring 
of antibiotics, the mutational alteration of the drug target to 
reduce the target's affinity for the drug or to reduce the 
permeability of the cell envelope. One example of the last one 
is the use of active transporters that belong to the ATP binding 
cassette (ABC) superfamily and utiUze the free energy of ATP 
hydrolysis to expel drugs from the cell (39, 47). Two unique 
tester sequences that could codify proteins involved in the 
antibiotic resistance of B. elkanii were found in this work. They 
were the sequence of clone 165 that was similar to the gene that 
codes for the peptidase M19, which has the abihty to cleave the 
P-lactam ring of penem and carbapenem antibiotics hke 
carbenicihn, and the sequence of clone 15 that was similar to 
the gene that codes for toxin secretion ABC transporter {hlyB 
like protein) of Bradyrhizobium sp. ORS278, which could be 
involved in the exporting of noxious molecules from the 
bacterial cell. 

The two remaining unique B. elkanii sequences identified 
by the RDA approach were both classified as genes that coding 
for conserved hypothetical proteins with no function predicted 
yet. One of them (clone 1170) was similar to a DNA region 
from Rhodopseudomonas palustris genome and the other 
(clone il57) had showed a similarity of 76% to a genomic 
DNA region of Nitrobacter winogradskyi. 

The possibility of RDA technique to identify unique 
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sequences had been confirmed in this work, since several DNA 
fragments unique to the tester DNA were obtained. These 
fragments included a large number of sequences that were not 
present in the genome of B. japonicum USDA 110. Some of 
the genes identified could be involved in mechanisms related to 
antibiotics resistance, motility, photosynthesis and transport, 
and could be the responsible for the phenotypic differences 
found between these two bacterial species. Further studies will 
be necessary to confirm the expression and functionalities of 
the predicted gene products. Nevertheless, the information 
provided will contribute to our knowledge concerning the 
molecular differences between these two species as well as the 
B. elkanii unique sequences identified could be used as DNA 
markers in order to investigate the rates of genetic material 
gain and loss in natural Bradyrhizobium strains. 
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